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Abstract 

An effective residual interaction for a meson-meson system is computed in lattice 
QCD. We describe the theoretical framework and present its application to the 1 = 2 
channel S-wave interaction of the vr-vr system. Scattering phase shifts are also computed 
and compared to experimental results. 



1 Introduction 

For many decades boson exchange models have served as a theoretical basis for mod- 
eling the strong interaction [Q. This class of models has reached a high level of so- 
phistication and has become a standard tool encompassing a large variety of strong 
interaction phenomena at low and intermediate energies 0. On the other hand it has 
also become increasingly clear that quantum chromodynamics (QCD) provides the 
microscopic underpinnings for hadronic interaction. 

The strong nuclear force operates in an energy range where the intrinsic hadronic 
ground states and a few excited states are important. At this low-to-intermediate energy 
scale the physics of QCD is only accessible through nonperturbative methods. While 
quark models have been used to capture part of the microscopic physics 0, there is 
no substitute for exploring hadronic interaction on a fundamental level. Currently the 
most appropriate way to do this is lattice field theory, supplemented by computational 
techniques. 

As of today lattice QCD is well developed and explores many different topics rang- 
ing from quark confinement to hadron masses. Hadronic interactions, however, still 
have not moved into the mainstream of activity. Early work of Richards and collabo- 
rators 0, Q addresses the role of quark exchange diagrams for meson-meson systems. 
Liischer's formula P], which exploits finite-size effects, allows the extraction of scatter- 
ing lengths from the lattice. It was applied to systems of pions and nucleons within 
lattice QCD. Beyond static properties of two-hadron systems theoretical concepts be- 
come more involved. Liischer has devised a method ^ for relating the excitation spec- 
trum of a two-body system in a finite box to scattering phase shifts. Elastic channel 
resonances are also within its scope [^. Phase shifts become available for a discrete 
set of momenta. Discrete interpolation between those can be done through the use of 
different-sized lattices. Within an 0(4) symmetric 0^ model extensive numerical tests 



T0[] have shown that the method works well. However, unlike with the previous model, 
stringent numerical standards are much harder to meet for applications with QCD, 
where gauge fields and composite asymptotic states complicate the issue. In Liischer's 
method the 'master' equation that relates the two-body energy spectrum to scatter- 
ing phases, unfortunately, has multiple solutions, the physical ones can be difficult to 
discern by solely numerical means. 

In the nonrelativistic limit interactions can be described by potentials. Ongoing 



lattice work by Green and coworkers [|T^ focuses on geometric configurations of quarks 
in hadron-hadron systems and their interaction energies. A rather large body of work 
exists |Tl|]-[|T^ with SU{2) as the gauge group, including the heavy-light meson-meson 



system studied by Steward and Koniuk ||T7|. 



In SU{3) a simple geometric model with static quarks tried by Rabitsch et al. 
T8|| provides some exploratory insight into the gluon-generated part of the interaction 



between two three-quark clusters. Extending this hne of work, studies of heavy-hght 
meson-meson systems were done |1^ where the heavy quarks are treated as static and 



the hght valence quarks are dynamic in the sense that their propagation is computed 
from the fermionic, staggered, lattice action. Although these studies are numerically 
difficult it appears that an attractive force between the partners is a common feature. 
However, in line with the exploratory flavor of these studies, the assignment of real- 
world quantum numbers to the hadronic systems was neglected. 

Two-hadron systems with one inflnitely heavy quark in each hadron have the ad- 
vantage that their relative distance r is a theoretically well-deflned notion. This opens 
the door to interactions of heavy-light systems. Still, realistic cases, involving B mesons 
for example, do require massive computing power. Most recently the UKQCD group 
has taken up simulations in this direction pO |. 



For light-light systems some work has been done in dimensions smaller than d+1 = 4 
and gauge groups simpler than SU{3). Probably the simplest nontrivial system in this 
class is a planar, d = 2 lattice theory with a U{1) gauge group. In flnite volume, and 
with properly chosen coupling, /3, the theory is conflning. It can be used as a 'toy' 
system to develop techniques for studying hadronic interaction. This was done within 



the staggered fermion scheme for a meson-meson system |21|. In the nonrelativistic 
approximation a potential can be computed and used in a Schrodinger equation to 
calculate scattering phase shifts. 

We will here follow essentially the 'blueprint' of ^^ applying it to a tt-tt system 
in the isospin 1 = 2 channel. Using an SU{3) gauge group in d = 3 dimensions and a 
highly improved action we attempt a realistic calculation of the residual tt-tt long-to- 
intermediate range interaction. Extrapolation to the chiral limit and a comparison of 
the corresponding scattering phase shifts to experimental results are performed. 



Preliminary stages of this work have been reported before [^, ^. A self-contained 
presentation of the formalism, in section |^, and a discussion and interpretation of the 
lattice results, in section ^, are the purpose of this paper. 

2 Formalism 

In the context of scattering within the framework of a quantum fleld theory the LSZ 
reduction formalism comes to mind. Putting the concept of 'in' and 'out' states, as 
t -^ ±oo in Minkowski time, to use in the Euclidean formulation proves, unfortu- 



nately, prone to theoretical difficulties. Maiani and Testa ||2^ studied a 3-point function 
{(t)qi(ti)(f)q^(t2)J{0)) involving a local source J and two fields (f) interpolating a pseu- 
doscalar particle. They showed that in the region -C ^2 ^ ^i? of Euclidean times, the 
3-point function depends only on the average (sum) of 'in' and 'out' matrix elements. 
Relative phases containing physical scattering information can thus not be extracted. 
This discouraging result, sometimes referred to as a no-go theorem, does not mean, how- 



ever, that it is in principle impossible to extract physical two-body information from 
the lattice. The no-go theorem can be avoided in special cases |2^, and also bypassed 



entirely by using a different method. A case in point is the aforementioned formal- 
ism by Liischer |p. It is based on a standard mass calculation involving the two-body 
system in a finite box, say with periodic boundary conditions. The discrete two-body 
energy spectrum obtained from a lattice simulation comprises the residual interaction 
provided the box is large enough so that the particles (hadrons) can separate from each 
other. In Liischer's original approach the spectrum of the interacting two-body system 
is in some sense 'compared' to that of free particles in the same box, their propaga- 
tion being represented by 'singular periodic solutions' of a Helmholtz equation. The 
method employed in the present work is similar in spirit to Liischer's to the extent that 
it is based on a comparison of two-body excitation spectra in a finite volume. We will, 
however, analyze the two-body information from the lattice in different ways. First, 
instead of using 'singular periodic solutions' to represent the noninteracting system, 
we construct a two-body system of noninteracting composite mesons within our lattice 
simulation. Second, we aim at a (nonrelativistic) potential between mesons as a means 
of interpolation between the lattice sites. 

In this section we describe the strategy for computing an effective residual interac- 
tion between composite hadrons from the lattice. The formalism will be set up for the 
example of two pseudoscalar pions, however, extensions to other hadron-hadron sys- 
tems should be straightforward: A set of interpolating fields for the two-hadron system 
is used to compute a time correlation matrix for the full (interacting) system. It has a 
time correlation matrix for the free (noninteracting) system as an additive term. Com- 
parison between those two correlators contains information about the effective residual 
interaction^. 

2.1 Meson-Meson Fields 

Starting with the one-meson field we choose 

Ut) = L-' E e^'^-' 4'2 (^, t) 75 4'2 ix,t), (1) 

X 

where L^ is the space lattice volume, and tp\^ {x,t) is a linear combination of quark 
fields '^pfA^J.{x, t) with flavor /, color A, and Dirac index /x, at lattice site x = (x, t). The 
field is local, however, the superscript ^^^ indicates that the field is 'smeared', iterated 
S times. The purpose of smearing p6|] is to enhance the amplitude for propagation of 
the hadron in the ground state. 

Products of those fields probe for multi-meson excitations. Our choice is 

%it) = 4>,pit) </)+p.(t) . (2) 



""^We closely follow [Ell in this section. 



These operators describe two-meson systems with total momentum P = and relative 
momenta p. They are designed to excite degrees of freedom of relative motion of the 
two composite mesons, in their respective intrinsic ground states. On a periodic lattice 
the momenta are 

2,71 -> -> Q 

p=—k with A; G R^ (3) 

2.2 Correlation Matrices 

Correlation matrices describe the propagation of the above fields in Euclidean time. 
Since we desire a comparison of the free (noninteracting) and the full (interacting) 
systems two types of correlation matrices are needed. 

The time-correlation matrix for the one-meson system (2-point function) is 

c^f^t^to) = {4{t)Mto)) - {4it))(Uto)) , (4) 

with ( ) denoting the gauge-configuration average. In the usual way it is readily worked 
out by means of Wick's theorem 

n n n 

. . . V^^^> (x) ^{^> (y) ... = .. . G^^} {x,y)... , (5) 

where the pair nn on the l.h.s. defines the partners of a contraction. We assume the 
quark propagator G^^^ of the smeared fields to be flavor-independent 

GfAf,,gBui^^ y) = ha GaIbu^x^ y) ■ (6) 

Since the smearing prescription used |2^ does not mix flavor indices, @ is a conse- 
quence of the corresponding relation of the unsmeared flelds and thus holds indepen- 
dently of 5*, see Appendix 0. We therefore will omit ^^^ from G to simplify notation. In 
(^ the separable term is zero, while there appear two contractions, n = 1,2. Making 
use of translational invariance one obtains 

with 



Cmt,t,)=6,^,^c,it,t,)e'P-^\ (7) 



Cp^(t,to) = L-3^e-*P^-^^(|G'A;.,B.(xt,foto)H • 



X 



Sums over repeated indices are understood. The source point xq = (xoto) is flxed and 
arbitrary. For later purposes we note that 

c_j?(t,to)=4(t,to). (9) 

The full meson-meson system propagates according to the 4-point function 

C^%to) = (^(t) %it,)) - (4(t))($,-.(to)) . (10) 

5 



Working out the contractions here leads to a more comphcated expression 



c^it,to} = ^-"EEEE 



^ip- {x2 -xi)+iq- {y-i - j/i ) 



^i ^2 3/1 2/2 
^*A2ti2,B2V2 (^2i^, y2h) GA2fi2,B2U2 {^2t, mto) X 

+ G\^^_^^^^^^{xit,y2to) GAi^ii,B2U2{^it,y2to) x 

GA2fi2,Biui (^2^, yito) GA2^l2,BlUl {X2t, yito) 
~ GA2fi2,Biuii^2t,yito) GA2fM2,B2V2{^2t,y2to) X 

Gauii,B2U2 (^l''-' ^2^o) Caiaji.Bii^i (a^l^^; ^l^o) 
-GA2^i2,B„yA^2t,yito) G*A^^^B^^^{x2t,y2to) X 

GAifii,B2V2{^l't,y2to) G\^^^^^_^^^{xit, yito) 



111 



Again, the separable term in (|1^) vanishes. Figure |T] shows the diagrammatic classifi- 



(xt) (x,t) (xt) (x,t) (xt) (x,t) (xt) (xt) 






(y,t,) (y,t,) (y,g (xt,) (y,t,) (y^g (y,t,) (yt,) 

A B C D 



Figure 1: Diagrammatic classification of the full correlation matrix C^^\t,to) according 

to fin-fO). 



(12) 



cation of ([n|), leading to four terms 

(j{4) ^ (j{4A) j^ (j{4B) _ (j{4C) _ (j{4D) 



It is useful to trace the pattern of contractions back to the 'meson content' of the 
correlators. Denoting the contractions as in (Bl) we write 



43 21 12 34 43 34 21 12 



21 43 12 34 43 34 21 12 



C(^B) 



^+p^-p^-g>+g) = (</'-p^+g^+p^-g 



U-,) ■ 



(13) 
(14) 



Each pair of equal numbers n = 1 ... 4 identifies contracted quark fields as they emerge 
from (|l|). Diagrams C and D exhibit valence quark exchange between the mesons 

23 41 12 34 

C^'""^ = (0l//>V-,>+,) (15) 

41 23 12 34 

C(^^) = (0 W_,-^+,) , (16) 

and thus must be considered as sources of effective residual interaction. 

We wish to extract from the full correlator C^^-* a part which describes two compos- 
ite 'lattice' mesons with their residual interaction 'switched off'. This part is entirely 
contained in diagrams A and B, but gluonic correlations between the mesons are also 
present because the gauge configuration average ( ) is taken over the product of all 
four fields. Gluonic correlations contribute to the effective residual interaction. It is 
thus necessary to isolate the uncorrelated part contained in C*^^"^-* + C^^^\ 

2.3 Free Meson-Meson Correlator 

Towards this end consider C^^"^-* in the form of (|T3|). In the spirit of (|^) we have 

^(4^) ~ {G*G G*G ) , (17) 

where the ~ indicates the Fourier sums, etc., which carry over from (|I|), see (pHI). The 
gauge configuration average in ([l^) may be analyzed in a systematic manner by means 



of cumulant expansion ||2^ . Taking advantage of {G) = we have 






4321 43 21 41 23 42 31 4321 

{G*G G*G ) = {G*G ){G*G ) + {G*G ){G*G ) + {G*G*){G G ) + {{G*G G*G)) . (18) 

The last term defines the cumulant. The first three terms on the right-hand side of 
(P^) are illustrated in Fig. B 



4321 fx \ h \ 4123 \ J \ 42 31 

<G*G><G*G> = t 4 T 4 <G*G><G*G> = Y i Y I <G*G*><GG> : 



Figure 2: Illustration of the cumulant expansion of C^^^' . The solid lines represent 
quark propagators, and the dashed hues gluonic correlations. 

The dashed lines indicate that the quark propagators are correlated through gluons. 
Evidently only the first one of the three separable terms in (|I8|) represents free, uncor- 
related, mesons. All other terms are sources of effective residual interaction between 

7 



the mesons. We therefore define 

-UA) 1 '' '\ '' 4 3 2 1 

C^ ^ = i^U^^) (0 V-.^) ~ iG*G ) {G*G ) . (19) 

A similar analysis of C^'^^^ leads to 

43 34 21 12 

c'^'''^ = {<pl,^-,^){<Plp<P-^)- (20) 

The sum of those is the free meson-meson time correlation matrix 

Cf^it^to) =C^^^\t,t,) + C%^\t,to) , (21) 

which is an additive part of the full 4-point correlator 

CW = C^') + Cj") . (22) 

The remainder Cj , defined through (|2^), comprises all sources of the residual inter- 
action, be it from gluonic correlations, or quark exchange (or quark-antiquark loops, if 
the simulation is unquenched). 

(4) 

The free correlator C describes two noninteracting identical lattice mesons. Their 
intrinsic structure is consistent with the dynamics determined by the lattice field model 

— (4) 

and its numerical implementation. Naturally, C m 
2-point correlator C'-'^\ Using (^ and (0),(^) gives 



-(4) 

and its numerical implementation. Naturally, C may be expressed in terms of the 



■^(i) _ ^(2) ^(2) ^(2) ^(2) ,^o\ 



pq ^p,q ^-p,-q ' ^-p,q^p,-q ■ 

Continuing with (|^ and (|^) this becomes 



Cf^it, to) = (<5p-,,- + 5_p-,,-) \c,it, to)f . (24) 



It is interesting to note that the property 



7^(4) _ 77(4) _ 77(4) _ 77(4) ...c-N 



which is evident from (pHf ), and also holds for C*-^-* through inspection of ([TTD , refiects 
Bose symmetry with respect to the composite mesons. Permutation of the mesons, one 
with momentum -|-p the other with momentum —p, results in the substitution p — > —p. 

2.4 Residual Interaction 

Eventually it would be desirable to extract scattering amplitudes, t-matrix elements, 
directly from a numerical simulation of lattice QCD. This would require a formulation of 
the LSZ prescription appropriate for the Euclidean lattice with the added complication 
that the asymptotic states are composites which themselves are determined by the 



discretized quantum field theory and its specific numerical implementation. We are 
presently not aware of such a formulation. 

For the time being we are content with defining an effective interaction from the 
computed lattice correlation matrices C^^-* and C along the lines of Appendix 0. 
We have in mind a description of the meson-meson system in terms of an effective 
Hamiltonian 

n = no + ni, (26) 

where Hq describes free propagation and Hi the residual interaction. We here adopt 
a definition of Hi suggested by an elementary Bose field theory subject to canoni- 
cal quantization. In such a theory it is straightforward to calculate the perturbative 
expansion of C^^\ as defined in ( p!0D together with (^, but using elementary Bose 
fields instead. This is outlined in Appendix pi The zero-order term of the expansion 

— (4) 

corresponds to C . Working up to second order we are lead to the following strategy: 
Define an effective correlation matrix 

C{t,to) = C^'\t,tof'^'c^'\t,to)C^'\t,tof'^' (27) 

= C^'\t,to)~'^' ci'\t,to) C^'\t,to)~"\ Bhhl , (28) 

where (^) was used for the alternative form.[| Then in terms of C we define 

t^^ dt ^ ' 

To obtain Hi numerical diagonalization of the effective correlator matrix C is called 
for. The eigenvalues of C are expected to behave exponentially with t, thus rendering 
Hi a time-independent matrix. 

In the numerical simulation it is desirable to utilize time slices with t as 'early' 
as possible in (^) since this gives smaller error bars. Towards this end the standard 
practice is to construct the interpolating fields from 'smeared' operators, see Appendix 
^. In this way the time correlation functions assume their asymptotic behaviour already 
at small t. 

Another technique towards the same end was described in [E3], and adopted in 



Ton foi' a lattice scattering problem. There, using the notation of |T^ for the moment, 
the generalized eigenvalue problem C(t)w'^ = Ai/(t, to)C(to)w'^ is considered. It leads 
to study the matrix V(t,to) = C~2(to)C{t)C~2(to). This is an interesting parallel to 
the form of (^7|). It is stated in |l^ that a reliable determination of the energy levels 
already for smaller values of t is achieved. 



2r^.-... . .- ..__._____,. ..__^ -(4) 



Dividing out square roots of correlators C for composite free mesons bears some similarity to 
'amputating' dressed external propagator lines. 



2.5 Lattice Symmetry 

From a numerical point of view computing C is greatly facilitated by utilizing symme- 
tries of the lattice. The action used with our L^ xT lattice is invariant under the group 
0(3, Z) of discrete transformations of the cubic sublattice. Using common nomencla- 
ture [^ the irreducible representations are T = Af,Af,E~^,T^,T^ with respective 
dimensionalities N-p = 1, 1, 2, 3, 3 . 

Given a fixed discrete lattice momentum p', see (^), the application of all group 
transformations g G 0(3, Z) generates a set of momenta p = Ogp' which all have the 
same length p = \p'\. These transformations define a representation of 0(3, Z), say its 
basis vectors are |p >, which is in general reducible. Let 

|P >= EE l(r,P)e X {T,p)e\p>, (30) 

r t 

with I (r, p)e >, e = 1 . . . A^r, being a set of basis vectors of the subspace that belongs to 
r. Systematic construction of those basis vectors is straightforward p9|. For example, 
choosing an on-axis momentum p' = {p, 0, 0), and F = Af , one obtains 

|(A+,0)1> = 1(0,0,0) > forp = (31) 

\{At,p)i> = ^E(l(±P'0.o)>+l(0'±p.o)>+l(o.o.±p)>) 

for p>0. (32) 

It is also useful to know that representations £ of the continuum group 0(3) char- 
acterized by their angular momentum £ = 0, 1 ... 4 have decompositions 



Q = 


= At 


1 = 


- Tr 


2 = 


= E-^®T+ 


3 = 


— 74-2 ^7 -t ]_ (i) -/ 2 


4 = 


= A+ © ^+ © T+ © T/ 



(33) 

if £ > 4 the last four lines apply cyclically. The above follows from letting 0(3, Z) 
operate on harmonic polynomials r^Yimi^, <P) which form a basis for L 

Since both the full and the free correlators O*^^-* and C , respectively, commute 
with all group operations g G 0{d^ Z) there exist reduced matrices, say C^^''^^ and 

— (4T) 

C ' , within each irreducible representation F such that 

< (F,p)e|OW(t,to)|(r',g)e' >= Srr'S.,Clf\t,to) (34) 

and similarly for C ' . It is obvious from (^) and (pSD that the reduced matrix 
elements of the free correlator have the form 

cfp{t,to) = 6pg |cf (t,to)f . (35) 

10 



The functions c^\t,to) are related to Cp'(t, to) by a p- and F-dependent factor. The 

— (4;r) 

above representation of C ' is of great advantage for numerical work because the 



inverse-square-root operation needed to compute the effective correlator (pTf ) is now 
trivial. In the sector F we simply have 

Cif HMo) = ^nT^S^^^^^^?^— • (36) 



-jn 



cip^\t,to)\\c'g'it,to) 

We envision numerical diagonalization, say 



Nr 

Ci,f\t,to) = E vPiP)>^'n\t,to)vP{q) . (37) 

n=l 



For asymptotic times the eigenvectors v^^ are time independent ||28[, whereas the 
eigenvalues behave exponentially 

Af (i^,^o) ~af exp(-u;f)(t-to)). (38) 

It is worth noting that diagonalization on each time slice, if numerically feasible, en- 
sures that excited states decouple from the ground state, and each other, by virtue of 
orthogonal eigenvectors. Thus, at least for a large enough size of C^^'^\ the t — > oo 
behavior avoids the two-body relative ground state. 



Extracting the effective residual interaction as defined in ( pQ]) now becomes a trivial 
matter. Provided that the v^^ are orthonormal we have 

Nr 

nfl = j:-iPip)^P-iPi^)- (39) 

n=l 

These are the desired matrix elements of the effective residual interaction in the sector 
F. In the basis \p> of lattice momenta we have 

<p|^/|g>=EE<pl(r,p)e>M!p^,<(r,g)e|g> . (40) 

r e 

How many irreducible representations F of the cubic lattice symmetry group actually 
are needed in the above series depends on the hadron-hadron system under consider- 
ation, the kinematical conditions, and on the partial waves that couple into the total 
spin J of the system. In this context the decompositions (Q) become relevant. 
In the numerical part of this work only the sector F = Af is considered. 

2.6 Nonrelativistic Potentials 

From a kinematical point of view a tt-tt system is ill suited for a description by non- 
relativistic potentials. Despite to-be-expected strong relativistic effects there are still 

11 



good reasons, however, to pursue this venue. First, we expect a potential to give insight 
into the nature of the residual interaction (attractive, repulsive, range, etc.). Second, 
a lattice simulation typically is performed in the realm of large pseudoscalar mass, 
certainly larger than the physical vr mass. In this region the potential picture still may 
have some validity. It is only through extrapolation that the relativistic region will be 
reached. Third, it seems that the shortage of lattice work in the area of hadron-hadron 
interaction, particularly for real- world systems, lends some justification to drastic ap- 
proximations, at this stage. 

The coordinate space matrix elements of the effective residual interaction are ob- 
tained from (|40|) by discrete lattice Fourier transformation. After redefining momentum 
variables we have 

< f 1 7^/ 1 s > = L-^ 5] 5] e'P- ^''-^e-''^- ^''^^ < p - q\ni\p + q > . (41 ) 

P Q 

It is useful to write the matrix element of Tij as a sum of < —qlTiil + q > and a 
remainder. A reference system where the relative momenta before and after a scattering 
event are ±q, respectively, is known as the Breit frame |^|. Then, the sum over pin (|4l|) 



gives rise to 6fs- and thus to a local operator V(f ) for the Breit-frames contributions, 
and a genuinely nonlocal operator VV(r*, s) for the remainder 

<f\Hi\s>=6^^gV{f) + W{f,s). (42) 

Effectively, V and W comprise all microscopic dynamical effects, ultimately by way 
of the lattice action, that possibly contribute to the residual interaction between the 
composite mesons. Loosely speaking, these include potential as well as kinetic (hopping) 
effects. In a quantum field theory this distinction hardly makes sense, however, there 
exists an interesting parallel in the many-body theory of nuclear reactions. There, 
techniques known as resonating group method (RGM) and the essentially equivalent 
generator coordinate method (GCM) lead to effective interactions between composite 
nuclei which are built from both kinetic and potential energies of some microscopic 
Hamiltonian. The brief digression of Section |2.6.1| will illuminate the analogy. 

Following common semantics we will refer to V and W as potentials. Their particular 
forms are 

V(f) = J2e-''''^-''<-q\ni\+q> (43) 



wif,s) = L-^EE 



gjp-(r-s)g-j(j-(r+s) 



P Q 



{<p-q\l-Li\p + q> - <-q\Hi\+q>) . (44) 

We will here limit ourselves to the relative S-wave local potential. Projecting (|43|) 
onto the partial wave ^ = 0, 

Vo(r) = ^y'rfa-V(r), (45) 
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implies a prescription for interpolation between the lattice sites. Whereas f in ( ^Bj ) is 
assumed discrete it has been reinterpreted as a continuous variable for the purpose of 
integration in (^Sj) . The effect of the angular integration on (^3|) is that only diagonal 
Af sector reduced matrix elements of the expansion (^) survive. We have 

Voir) = Y.U2qr)n['^l\ (46) 

Q 

where jo{x) = sin(x)/x. Note that the sum in (|46|) runs over discrete lattice momenta. 

2.6.1 Succinct RGM 

We divert to specify the analogies mentioned above. The resonating group method 
(RGM), and it's cousin the generator coordinate method (GCM), are specifically ap- 
plicable to a system of two composite nuclei. For the many-body wave function of 
two nuclei, say 1 and 2, the RGM ansatz is $ = A {(l>i{Ci)'p2iC2)gix)) , where ^1,2 are 
intrinsic coordinates and x is the relative coordinate between the clusters. The oper- 
ator A antisymmetrizes with respect to nucleon exchange. The microscopic Hamilton 
operator H, with AH = HA, can be written a.s H = Hi + H2 + T12 + Vu where 
Hi 2 refer to the intrinsic degrees of freedom and T12, V^i2 are the relative kinetic and 
potential energy operators. While 0i^2 are input, for example shell-model wave func- 
tions, the relative wave function g{x) is subject to solution of the RGM equation 
{(l)i4)2Sff\{H — E)A\(f)i(p2g) = 0, all r, where Sff{x) = S^{x — r). In terms of the in- 
tegral kernels n{f,r') = (0i02'^r|-^|0i02'^r') and h{f,r') = (0i02'^r|-f^-^|0i02'^r') the 
RGM equation reads {h — En)g = 0. Assuming that the so-called overlap kernel n has 
no zero eigenvalues (redundant states) it can be cast into the form of a Schrodinger 
equation. Defining 

hnGM=n~'/'hn-'/\ (47) 

we have {hjicM — E)f = with / = n^/'^g. Here hjicM = h-Q + hj is an effective 
Hamiltonian for the two- nucleus system. Its free part, ho, stems from splitting off 
terms with no rest antisymmetrization, direct terms, in this case the first terms in 
A = 1 - {1 - A) and HA = {Hi + H2 + T^) + V12 + H{A - 1) in the definitions of 
n(r, r ') and /i(r, r '), respectively, and also in n"^'"^ = 1 + (n^^'"^ — 1). Direct terms lead 
to local operators on the wave function /(r). In this way (Jirgm — E)f = assumes 
the form 



h' 



2m 



12 



\-f{r) + VRGM{f^ f{r) + 1 d'r'WRGM{r,r') fin = {E-Ei-E2)f{r). (48) 



Here 77112 is the reduced mass, the kinetic operator is identical with T12, and Ei^2 are 
the internal cluster energies. Microscopic dynamics with kinetic and potential energies, 
contained in H, are cast into the effective potentials Vrgm, Wrgm- 
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Analogies with ( P^D and (^) suggest themselves. Since the spectrum of Ho is akin 
to 2\/vn,^ + p^ = 2?7i + p^ /m + . . ., similarities exist between Ei + E2 — 2^— A^. and 
Ho, and between Vrgm, Wrgm and V, W or 7i/. Loosely speaking 2m plus the relative 
kinetic energy of our two-pion system drops out in the ratio of (|27D. 

3 Computational Issues, Analysis and Results 

The numerical implementation poses some extraordinary challenges. Compared to typ- 
ical hadronic mass scales, energy shifts due to residual interaction can be quite small. 
On the other hand, the total energies of a two-hadron system at larger relative momenta 
tend to be big. Thus extracting energy shifts from ratios, see (|36|), of steeply dropping 
time correlation functions can push the numerical quality of the lattice simulation to 
its limits. 

3.1 Lattice Implementation 

The computation has been done on an L^ x T = 9^ x 13 lattice. We have used a tree- 
level Oi^a"^) tadpole improved next-nearest-neighbor action, as discussed in [^, ^, in 
a quenched simulation. The number of gauge configurations was Nu = 208 separated 
by at least 1024 update steps. The hadron spectrum of this action has been studied in 
|33[] , and also in [0], on coarse lattices. Quite good agreement with the experimental 



masses was observed in both cases. In the conventions of [33] the gauge field coupling 



parameter was /5 = 6.2. This corresponds to a lattice constant of 

a ^ 0.4 fm or a"^ ~ 0.5 GeV (49) 

as determined from the string tension^. This lattice is rather coarse, on the other hand 
a large lattice volume (here La ~ 3.6fm) is required to accommodate two hadrons and 
to allow for sufficiently large separations at which their residual interaction becomes 
negligible. 



The fermion part of the improved action |33[ contains the Wilson hopping parameter 



K. Quark propagators were obtained for the six values 

K-^ = 5.720, 5.804, 5.888, 5.972, 6.056, 6.140 , (50) 

using a multiple mass algorithm [^]. The critical value for k,~^ is about 5.5. The 



correlator matrix (|TT|) requires quark propagator matrix elements G{xt, yto) between 
arbitrary lattice sites, except that the initial time slice is fixed. Computing all of those 



is not feasible. A random source technique very similar to the one described in [36 



^ Unless otherwise noted all lengths and energies are given in units of a and a ^, respectively, 
throughout the text and the figures. 
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was employed, see Appendix ^. The number of complex Gaussian random sources was 
Nji = 8 for each color-Dirac source point, thus we effectively have 8 x 3 x 4 = 96 
random sources per gauge configuration available to estimate G. 

The fermion fields were subjected to 'Gaussian smearing'. We have used the smear- 
ing strength parameter a = 2, in the notation of |2^, and 5 = 4 iterations for smearing 



at the sink. This parameter set is similar to the one used in pSl with the same type 



action. In addition 'APE fuzzing' ^7j , which is the analogue of smearing for link vari- 
ables, was done on each gauge configuration prior to smearing, also with a = 2 and 
S = A. 



3.2 Correlator Analysis 

The centerpiece of the numerical effort is the effective correlator matrix ([361). We have 
computed 5x5 A^-sector matrices C^^'^'^ \t, to) and C ' ^ (t, to) for on-axis momenta 



277, 



Q 



-k with A; = 0,1,2,3,4. 



(51) 



The source time slice is at to = 3. It turns out that the C^^''^'^ ^(t, to) are dominated by 
their diagonal elements. The off-diagonal elements are small and have somewhat large 
statistical errors, mostly extending across zero. Hence it is reasonable to adopt the 
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Figure 3: Free and full correlator functions (C) and their ratio (R), for k, ^ = 5.720 
and momentum k = 1, showing a repulsive mode. 
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Figure 4: Free and full correlator functions (C) and their ratio (R), for n ^ 
and momentum k = 3, showing an attractive mode. 



5.720 



diagonal approximation, replacing 



Thus the effective correlator becomes 



cr^^Mo) 



^(4;A+) 

Lyqq 



it,to) 



C 



<4;A+) 



qq 



it,to) 



(52) 



(53) 



with off-diagonal elements set to zero. 

In Figs. P and § we show two examples of free and full correlator functions and 
their ratios, as they appear in (^3]), both for the smallest inverse hopping parameter, 
K,~^ = 5.720. The examples are for momenta k = 1 and k = 3, respectively, see (|5lD . 
For the repulsive level the mass shift is quite small, which translates into a very noisy 
signal for the ratio. For the attractive level the mass shift is rather large, however, the 
fall-off of the free correlator function is very steep, which has the effect that only very 
few time slices are available for analysis. The straight lines in Figs. |^ and | come from 
linear fits to the logarithms of the correlator functions. Only data points marked by 
filled plot symbols were used in the fits. 



The lattice action used here, see Sect. |3]T|, exhibits ghosts. These are unphysical 
branches in the lattice-quark dispersion relation and are indigenous to highly improved 
actions pSl. The presence of ghosts can contaminate the signal on early time slices. 
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In our case this contamination is clearly discernible only at t = 3. The corresponding 
data points of the correlator functions have been ignored in our analysis. 

Quantum fluctuations manifest themselves in an overall renormalization factor that 
relates the free and full correlators at t = to 

C!t^^"\to,to)=Z'cff\to,to)- (54) 

The value obtained from the fits is Z"^ = 1.97(3), or Z ~ 1.18. It is the same for all 

K, a consequence of the quenched approximation. We have plotted Z~'^Cqq ^ (t,to) in 
Figs. ^ and ^ to improve readability. Note that Z^ has no influence on Tij. 



3.3 Residual Interaction 

In Fig. |] we show an example of the energy level shifts of the meson-meson system due 



W 3 



5.888 



free interactinc 



Figure 5: Energy level shifts (hyperfine splittings) due to the residual interaction for 
the momenta k = 0,1, 2, 3, 4 at a fixed k~^. 

to the residual interaction, for a selected k. Low-momentum (long-distance) modes 
are weakly repulsive, whereas high-momentum (short-distance) modes are strongly 
attractive. Errors stem from a bootstrap analysis. 

Analyzing the effective correlator (|53|) along the lines (37)- (39) yields the matrix 

elements Ti-j^g . These were used to compute S-wave local potentials Vo(r) according 
to (^61). For the current lattice the sum over on-axis momenta (PH), 



feniax 

q k=0 



(55) 
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truncates at fcmax = 4. We try to obtain some feel for systematic errors caused by 
the momentum cut-off /Cmax- On a coarse lattice one should expect those to be large. 
Thus, in Fig. ^ we show a family of potentials, all for k^^ = 5.888, which correspond 



> 



-2 - 



-4 




Figure 6: Fourier-Bessel representation (^) of the S-wave local potential, at k ^ = 
5.888. Various momentum truncations, see (|55|) , are shown. 

to various upper limits fcmax- 

The matrix elements of the residual interaction (|39|) depend on the pseudoscalar 



m 



mass m through the hopping parameter k. Figure shows Ti-j ^^ as a function of 
for momenta A; = ... 4. The solid lines are three-parameter fits using 

hq{x) = hq + h'qX + h'^x^'"^ with x = rn? , (56) 

which is motivated by chiral perturbation theory [^. Extrapolation to the chiral limit 

^SlP = h,{x ^ 0) (57) 

then yields matrix elements of Tii which describe the residual interaction of the physical 
vr-TT system in the 1 = 2 channel. We will exclusively refer to those in the subsequent 
discussion. 




Figure 7: Chiral extrapolations Ti.j gg (m? -^ 0) for momenta fc = . . .4, see (^T]). In 
physical units the six masses m range from 1.2 to 0.72 GeV, the filled plot symbols 
refer to 171^^=0.14: GeV. 



3.4 Parameterization 

In Fig. H the solid lines correspond to the extrapolated potentials for k^a.^ = 3 and 4, 
respectively. The oscillations caused by the representation as a truncated Fourier-Bessel 
series complicate their interpretation. 

To get a clearer picture, consider the parametric function 

1 - a2^"' 



\/(°)(r)=«i 



1 + a3r"5+^ exp(a4r) 



+ ao where r G [0, cxd) 



(5^ 



It is flexible enough to represent attraction and/or repulsion at certain ranges, and has 
the feature that the Yukawa form is approached asymptotically. 



y(")| 



r 



aia2 e 



-04r 



"3 r 
Now define a 'latticized' version of (|58|) via 



+ ao for r — > oo . 



Vrir) = E Jo{2qkr) l^f ^ with q, = —k 



27r 



(59) 



(60) 



fc=0 



.(a) 



The V^ are expansion coefficients which we have determined from calculating the 
matrix Bik = Jo(2gfcn), with support on tq = 0, n = j(k~lTi)^ 1 = 1... k^e.^, and 
then applying its inverse to (^) as 



vt^ = T.B^'y^"^^ 



ri) . 



(61) 
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Figure 8: Local potentials for momentum truncations A;max = 3 and /Cmax = 4. Shown 
are the raw lattice results in the m'^ -^ limit (extrapolated, solid line), the fits with 
the Fourier-Bessel trial potential V^" (r) (fitted, solid hne), and the corresponding 
r G [0,00) parameterization V^°'\r) (parameterized, dashed line). 



K") 



In this way the span of V£ (r) is a subspace of the space of functions Vo(r), see (^) 



and (|55|), accessible through the lattice simulation. Now make fits to lattice potentials 
minimizing 



L/2 



dr 



Vn(r) 



Vt\r] 



mm a 



(62) 



We have found that 05 = 2 gives an almost perfect match and thus held this parameter 
fixed. Also, Oq deviates very little from which indicates that the interaction 'around 
the world', across La = 3.6fm, is negligible. The results of the fits are the curves in Fig. |] 
marked 'fitted'. In fact the 'fitted' and the 'extrapolated' curves are indistinguishable 
within the line thickness. (This changes if a^ 7^ 2.) The curves marked 'parameterized' 
show the corresponding r G [0, cxo) parameterization V^"\r) as defined in (|58|). 



3.5 Physical Potential 

It is tempting to use the resulting fit parameters, shown in Fig. ||, to make contact 
with the boson-exchange picture of strong interactions, see Fig. p. This can be done 
identifying the Yukawa asymptotics (|59|) of V^"(r) with 



m^ 



a\ a 



1.4-0.9 GeV 



(63) 
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|! = -^1^^17.4-2.5, (64) 

An as 

see (1^). The left and right numbers above relate to k^ax = 3 and k^ax = 4, respectively, 
indicating large systematic errors. Both the mass m^ of the exchanged particle and the 
vertex coupling, g ~ 14.8-5.6, are typical for a hadronic system. For example g is 
around 13.45 for the vr-N vertex that enters the N-N interaction. Statistical errors are 
also large for these quantities (see Fig. p!T|). 





Figure 9: Boson exchange diagram related to (^) and (|63D , (|6^) 



The above discussion addresses the long-range physics of the system. Figure |T0 
shows the potentials in the region beyond the lattice cut-off, r> a, in physical units. 
Both truncations at /^max = 3,4 are in agreement and indicate a repulsive residual 



interaction. In Fig. |TI] we show the physical potentials for truncations fcmax = 3, 4, 
respectively, together with a family of potentials computed from 8 bootstrap samples. 
The spread of the dashed lines in Fig. |Tl| is indicative of the statistical fluctuations. 
Compared to the fcmax = 3 case the fluctuations of the fcmax = 4 potential are visibly 
larger. A generic feature of two-hadron systems comes to mind as a possible reason. 
Correlator matrix elements which involve large momenta describe very massive (two- 
body) states and consequently drop very steeply with t. On small isotropic lattices it is 
very hard to deal with the resulting deterioration of the lattice signal since data from 
only very few time slices are usable. From a numerical point of view we consider our 
results using the truncation fcmax = 3 the most reliable. Nevertheless, we observe that 
both momentum cut-off values lead to qualitatively consistent potentials in the long 
and intermediate range region, r > 0.4fm. 

The steep drop of the potentials from their values at r = a to deeply negative values 
at r = 0, see Fig. H, does not influence low-energy physics. It is an indication, however, 
that the physics in the system at r = is special. A possible explanation is suggested by 
the SU{3) color content of the two-body system. Using standard nomenclature [^, ^ 
we note that for the color structure of the one-meson interpolating field (|l|) only the 
singlet from 3®3 = 8©lis used. Thus the color-source structure of the two-meson 
interpolating field (0) is that of an overall singlet, 1 ® 1 = 1. A singlet is also contained 
in the decomposition of the product of two color octets 808 = 27©10©8©8©T0©1. 

21 



> 



0.6 



0.5 



0.4 



Q) 0.3 

o 



0.2 



0.1 - 



0.0 



-0.1 















- 


k = 

max 


4 

/ \ 

/ x 

\ 

\ 






- 


— 


k = 

max 


3 \ 


V 




— 


- 






\ \ 




- 


- 






\ 


\ 


— 


- 


k = 

max 


2 ,■■ 




^^.^^^ 


- 


- 










- 



0.0 



0.4 0.. 

r [fm] 



1.2 



Figure 10: Local potentials for k^ 



2, 3, 4 in physical units for distances r < 0.4fni. 



Naive gluon exchange in 8® 8 is attractive [0. However, in the confinement phase the 
propagation of the system into an 8 (S) 8 color configuration is dynamically suppressed 
with increasing distance r. This means that we should expect the interaction energy 
to be more attractive at r = as opposed to all other r = a, 2a, 3a . . .. The situation 
is of course more complicated because of dynamical effects from the spin (Dirac) and 
the flavor degrees of freedom. 

In simulations of heavy-light meson-meson systems performed in coordinate space 
19] the above mechanism is seen directly. There, the r = case can be easily isolated. 
In our Fourier-Bessel analysis, on the other hand, strong attraction at r = can bias 
the result at r ~ a via oscillations of the basis functions. The danger of course is that 
oscillatory features are misinterpreted as repulsion. We have tried to minimize this 
by employing the parameterization (pSD which would easily be capable of revealing 
attraction, say at r = a, but in fact gives a repulsive interaction in the intermediate 
and long distance region, r> a, as an answer. 
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Figure 11: Local potentials for momentum truncations fcmax = 3 and A^max = 4 (solid 
lines) and families of bootstrap samples (dashed lines) to estimate the overall statistical 



errors. 



3.6 Scattering Phase Shifts 

It is obvious that relativistic effects in the tt-tt system are very large in the experi- 
mentally relevant kinematic region. They are also inherent in the lattice simulation. 
Nevertheless it is interesting to calculate scattering phase shifts as they arise from 
the computed potential. This excercise is much in the spirit of section 2.6.1 , where it 
would be based on solving (^) for scattering states. We have used a Volterra inte- 
gral equation of standard potential scattering theory |^| with V^^^r) and obtained 



the scattering phases from Jost functions. The results are shown in Fig. |T^, where 
the two frames distinguish the truncations /cmax = 3,4 respectively. The two upper 
solid curves correspond to the physical reduced mass m7r/2 of the tt-tt system, with 
TT^TT = 0.28a^^ = 0.14GeV. In order to estimate the magnitude of the error due to the 
non-relativistic potential scattering theory we have also used the relativistic dispersion 
relation mj^^p) = J ml. + p^. The two lower solid curves in Fig. |l^ show the results. 

Estimation of the statistical errors on the phase shifts is not a straightforward task 
and requires some discretion. One option is to select pairs of extremal potentials from 
the bootstrap analysis of Fig. |n| and compute phase shifts for those. This was done 
separately for fcmax = 3 and 4. We have chosen the potentials which stand out in Fig. ^ 
as the curves with largest and smallest values, respectively, at around r = 0.6fm for 
^max = 3 and around r = 0.8fm for fcmax = 4. The resulting phase shifts are shown in 



Fig. 12 as the boundaries of the dotted regions. Not surprisingly, as discussed above, the 
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Figure 12: Scattering phase shifts 6^^^ calculated from the tt-tt potential of the lattice 
QCD simulation (thick lines). Influence of momentum cut-off k^a.x = 3 (left) and 
^max = 4 (right) is shown. Results using the classical and relativistic dispersion relation 
are distinguished by m,r and m.„{p), respectively. Errors are represented by the dotted 
regions. Their boundaries (dashed lines) correspond to the phase shifts calculated with 
extremal (bootstrap) potentials of Fig. |Tl|. The experimental data are from |^ . 



error spread is much larger for k^^x = 4. Results for both truncations are, nevertheless, 
consistent. 

In the Volterra equation the parameterized potentials V^^^r) were used for the 
entire interval r G [0,cxd). Their deep attractive trough within < r < a is probed 
by large momenta and is responsible for the scattering phase shifts turning upward 
beyond Mt^t^ > 0.8 — l.OGeV. This region is above the lattice cut-off. For smaller Mt^tt 
the intermediate and long distance region of the potential, essentially as displayed 
in Fig. |TU|, is probed. There, it appears that our results are in qualitative agreement 
with experimental findings. The data points shown as triangles in Fig. [l^ are from an 
analysis of the CERN-Munich vr+p — >■ vr+vr+n experiment by Hoogland et al. |^ . 
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4 Conclusions and Outlook 

During the last two decades QCD has emerged as the underlying theory of strong 
interactions. It is an important problem to put QCD to the test of explaining the inter- 
actions between hadrons. The nonperturbative nature of QCD in the energy domain of 
nuclear physics makes lattice field theory the predestined technique to deal with this 
question. 

We have made an effort in this work to develop and apply lattice QCD techniques 
to the problem of hadron-hadron interaction. Within the example of a meson-meson 
system we have outlined a strategy usable to extract a residual interaction from a lattice 
simulation, which also may point a way to application to other hadronic systems. 

On the practical side we have done a simulation of the 7r+-7r+ system and extracted 
a local potential for the S-wave interaction. A coarse and large- volume lattice to make 
space for two hadrons, the extraction of small shifts of energy levels, and steeply de- 
clining time correlation functions indigenous to two-hadron states all conspire to make 
the numerical work rather challenging. Nevertheless, our results indicate a repulsive 
interaction at long and intermediate ranges, down to r ^ 0.4fm, monotonously rising 
to ~ 0.4GeV at that distance. The corresponding scattering phase shifts are subject to 
large relativistic corrections in the kinematical region of interest, however, within those 
limitations, they compare favorably with experimental results. The parameterization 
of our lattice results, matching a Yukawa potential in the asymptotic region, allows us 
to make contact with traditional boson-exchange models. In this picture we find that 
the mass of an exchanged particle and its vertex coupling g are in the region of IGeV 
and 10, respectively, which are both reasonable for those quantities. 

It appears that lattice techniques as tried in the present work are feasible for tackling 
hadron-hadron interactions. Subsequent work should make use of anisotropic lattices 
to ease the problems related to steeply dropping correlation functions. It would also 
be highly desirable to go beyond the 'potential' picture and be able to extract the 
scattering amplitude (t-matrix) in a more direct fashion numerically from lattice QCD. 
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A Random Sources 



We here give details of the random source technique used to compute quark prop- 
agator matrix elements and of the smearing procedure. 
Consider the linear equation 

EE G-/^%ixx,,yy,)xi^;'''''''^<\yy,)=6^A'S,,,R^^'^'^^ (65) 



J/S/4 



Bu 



where G~^^^^ is the (known) fermion matrix for flavor /, and R are complex Gaussian 
random vectors of length L^ that live on the space sites x of the lattice. The meaning of 
the indices are A,B . . . = 1, 2, 3 color, /x, z/ . . . = 1, 2, 3, 4 Dirac, x,y... space (rf = 3), 
0:4,1/4 time, and r = 1 . . . Nn labels the random sources R for each source point. A 
prime ' denotes a source point. There is some freedom in choosing the latter. In (|65|) 
the sources are nonzero on one time slice x'^ = to only. A new source is chosen for 
each color, Dirac and time index. The same set of sources is used for different flavors 



(meaning n) in order to take advantage of a multiple-mass solver p5 



A version of a random source technique, called 'maximal variance reduction' |44 



IS 

currently in use for light-quark propagators in heavy- light systems |2y] . There, sources 
are spread across the entire lattice, including all time slices. Since the variance of the 
source typically is of order one, exponentially decaying time correlation functions will 
quickly be engulfed in noise. To alleviate this problem a subdivision of the lattice into 
disjoint regions, say < t < T/2 and T/2 < t < T, is made and noise reduction is 



achieved for propagators connecting those regions, see [^ . On the other hand random 
sources oc 6tto, as in (|65|), avoid the above problem from the outset for lack of noise at 
the sink. They have proven sufficient in ||21| so we continue to employ them here. 
The random sources R are normalized according to 

^ R^^'^'^<\x)R^^''''^y'^>{y) = 6a'b' V^' ^'^-^yl • (66) 

(r) 

Here J2(r) denotes the random-source average, which we approximate numerically as 



1 Nr 
(r) R r=\ 



Employing the solution vectors X of (|65| ) an estimator for the propagator matrix ele- 
ments then is 

G%\A^{yy,,xx,) = E^f£^""^^(m)i?(^^^^^)*(f) • (68) 

(r) 
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Operator smearing PBI is defined through 



ri°n^^) = ^a{x) ,lj^:\xt) =Y.Y.KAB{x,y) ij'^r'Hm , (69) 

B y 

with s G N, and the matrix 

3 

Kab{x, y) = Sab S^^g + a 51 \Um,AB{xt)5^^g-rh + Ul,ABiy^)^x,y+m ■ (70) 

m=l 

The real number a and the maximum value S for s = . . . S are parameters. We have 
used fuzzy link variables U G SU{3) in (ffOl). Due to the linearity of (|65D the above 
iterative prescription translates directly to the random source and solution vectors, 



4°»^'^"^*')(i") = 


6c'B'R^'''''"''''\z') (no sum over B') 


(71) 


4f(^'^"^*'ni") = 


> ,> ,^C7'B"(2: ,2/ )R]^„ '{y) 
y' B" 


(72) 


xS'^y-'^''"''\zt) 


^(f;B'u"rt')(^. 


(73) 


xijW;B'."rt')^^^^ 


X A 


(74) 



Finally, replacing R -^ R^^^ and X -^ X^^^ in (^) yields the propagator G^^^ for 
smeared fermion fields as used in (^ and @. 

B Elementary Bose Field 

Let 0(a;) be an elementary Bose field defined on the sites x = (x, t) of the lattice. 
It is understood that is subject to canonical quantization. Let £ = £o + £/ be a 
Lagrangian such that Co{(p, d(p) is the free part and £/ = Ci{(f>) is a (small) interaction. 
Thus C gives rise to a Hamiltonian H = Hq + Hj with according interpretation of 
its terms. In analogy to (H) and (|) further define 0p(t) = L^^ Y,x e^^'^(f){x, t) and 
$p-(t) = (p-pit) (t)+p{t) . We have in mind a perturbative calculation of the correlation 
matrix ^^ 

^?i(i,to) = (0|<|.J-.(t)<l,-(to)|0) (75) 

assuming a nondegenerate vacuum state |0) with ifolO) = 0. Standard perturbation 
theory gives rise to the time evolution operator 

oo l_l\N rt rt ^ 

Uit, to) = E ^iJr dti--- dtN T[///(ti) . . . Hi{tN)] , (76) 

^r_n A*! J ta J to 



N=0 
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where Hf(t) = e^°^^ *"-' Hje '^"*^* *"^ refers to the interaction picture. Working out ( [T5| ) 
with (|76|) induces a perturbative expansion of the correlator 



(4;iV), 



C'pqit^to) — Z2 ^p<f (^'^0 



(77) 



N=0 



Exphcit forms of the A^ = and the A^ = 1 terms are conveniently expressed in 
terms of wave functions ip^^J (p) defined through 



c^^XKp) = {ni^\K<to)\Or 



(7J 



where Inu) is a complete set of eigenstates of Hq with Ho\niy) = WJ['^\nu), and c^} 
are appropriate normalization factors chosen such that the ip!;^j} are orthonormal. The 
corresponding correlator terms are 



C. 



(4;Ar=0) 
pq 



{t,tc 



E 



(0)|2^-H.r(*-to)^(Oj(^-.)^(0).(^-.) 



(79) 



C^^''(t,t, 



pq 



{nu\Hj\mfi)c^^}*c^^lexp 






sinh 



[t — tQ)6nm + 



H-r-tyj") 



2 

{t - to) 



In *^ m 



(t- 


-to) 






1 


1- 


"nm 


) 



(80) 



Without loss of generality the normalization constants c^°j may be chosen real and 
positive. We now observe that the two normalization factors and the exponential in 
(pop may be removed by multiplying C(^'^='^) from both sides with the inverse square 
root of (7(^'^=°). Hence the matrix elements of 



-1/2 



C^^-'^='\t,to) = C'(^^^=°)(t,to)^ ' C'('^^=^n^,^o)C'(^^^=°)(t,to) 



-1/2 



(81) 



in the basis ip!;^J{p) are products of {niy\Hj\miJ,) and the expression inside {. . .} of (pO]). 
The t derivative of the latter is equal to one at t = to- Thus we have 



dC, 



{4;Af=l) 
pq 



it,to) 



dt 



-Y.Y. ^^nKp){nv\H,\m^^)^'^^>{q) = Hj 



(0)* 



(821 



t=to 



nv m/i 



which holds independently of the basis. Finally, we may replace C'*^^'^"^-' in (^)- (|8^) 
with the full correlation matrix C^'^^ since this will only introduce an A^ = 2 error. This 
leads us to define an 'effective correlator' 



C(^)(t,to) = C'(^^^=°)(t,to)~'/'c'(^)(t,to)C'(^^^=°)(t,to) 
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-1/2 



^3) 



We have shown at this point that C'-^^ has the power series expansion C*^^)(t,to) = 
1 + Hi{t — to) + 0{h'^) where the second-order remainder must depend on the product 
h = Hi{t — to) for dimensional reasons. In the next order N = 2 a. calculation similar 
to the above reveals that 

C^'\t,to) = 1 + Hiit -to) + l {Hi{t - to))' + Oih') . (84) 

Whether or not these are the initial terms of a converging power series expansion for 
C*^^^(t,to) will depend on the actual Hj. In case the series converges 'everywhere' the 
expansion will define a correlator even for large t — to- Limiting ourselves to order N = 2 
perturbation theory (|84|) is the same as 

C(4)(t,to) =e-^^(*-*»). (85) 

The utility of these results in the framework of a lattice simulation lies in the 
analogy that can be drawn between (7(^'^=°) and the free correlator C , and between 
C^'^^ and the full correlator C*-^-'. The analogue of (|85|) may then be considered as a 
definition of an effective interaction. 
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